Volume  99,  Number  3,  May-June  1994 

Journal  of  Research  of  the  National  Institute  of  Standards  and  Technology 

[J.  Res.  Natl.  Inst.  Stand.  Technol.  99,  247  (1994)] 

Three-Axis  Coil  Probe  Dimensions  and 
Uncertainties  During  Measurement  of 
Magnetic  Fields  from  Appliances 


Volume  99 


Number  3 


May-June  1994 


Martin  Misakian  and  Charles 
Fenimore 

National  Institute  of  Standards 
and  Technology, 

Gaithersburg,  MD  20899-0001 


Comparisons  are  made  between  the  av- 
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results,  which  are  determined  assuming 
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investigation  extends  an  earlier  treat- 


ment of  the  problem,  which  did  not 
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clusions reached  previously. 

Key  words:  appliance;  coil  probe;  mag- 
netic field;  measurement;  measurement 
uncertainty;  power  frequency. 

Accepted:  April  13,  1994 


1.  Introduction 

This  paper  reconsiders  a problem  related  to  the 
measurement  of  power  frequency  magnetic  fields 
from  electrical  appliances  using  three-axis  circular 
coil  probes.  Specifically,  it  reexamines  the  differ- 
ences between  the  average  magnetic  flux  density  as 
determined  using  a magnetic  field  meter  with  a 
three-axis  circular  coil  probe  and  the  magnetic  flux 
density  at  the  center  of  the  probe,  5o,  assuming  the 
field  is  produced  by  a small  loop  of  alternating  cur- 
rent, i.e.,  a magnetic  dipole.  The  “average”  arises 
as  a consequence  of  the  averaging  effects  of  the 
coil  probes  over  their  cross  sectional  areas  when 
placed  in  a nonuniform  magnetic  field.  The  differ- 
ences between  the  average  magnetic  field  and  iJo 
can  be  regarded  as  measurement  errors  because 
the  center  of  the  probe  is  normally  considered  the 
measurement  location.  The  magnetic  dipole  field  is 
chosen  as  the  relevant  field  because  its  geometry 
provides  a good  approximation  of  the  magnetic 
field  produced  by  many  electrical  appliances  [1]. 


The  average  magnetic  flux  density  measured  by 
a three-axis  magnetic  field  meter,  5av3,  is  also  re- 
ferred to  as  the  resultant  magnetic  field  and  is  de- 
fined as  [2] 


= + (1) 

where  Bi,  Bj,  and  B3  are  average  root-mean-square 
(rms)  magnetic  field  components  determined  by 
each  of  three  orthogonally  oriented  coil  probes. 

Differences  between  5av3  and  Bn  are  calculated 
as  a function  of  r/a  where  r is  the  distance  between 
the  magnetic  dipole  and  the  center  of  the  probe, 
and  a is  the  radius  of  the  three-axis  probe.  In  addi- 
tion, differences  between  5av3  and  Bn  are  examined 
for  different  orientations  of  the  magnetic  dipole 
and  rotations  of  the  three-axis  probe.  Because  the 
relative  orientation  of  the  dipole  and  three-axis 
probe  is  not  known  during  most  measurement  situ- 
ations, there  is  a distribution  of  possible  differ- 
ences between  5av3  and  Bn,  and  these  differences 
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collectively  represent  a source  of  measurement  un- 
certainty for  a given  rja . What  will  be  of  interest  in 
this  paper  is  the  largest  difference  that  occurs 
between  5av3  and  Bo  as  a function  of  rja  (for  all 
possible  orientations  of  the  dipole).  This  largest 
difference  is  designated  AZ^maxs- 
This  investigation  extends  an  earlier  treatment 
of  the  problem  which  considered  different  orienta- 
tions of  the  dipole,  but  not  all  possible  orientations 
of  the  three-axis  probe  [3].  The  maximum  differ- 
ence between  5av3  and  Bo,  A5max3,  is  found  by  a nu- 
merical search  during  which  J3av3  is  determined  by 
numerical  integration.  The  major  advance  over  the 
earlier  study  is  the  development  of  an  expression 
giving  the  average  magnetic  flux  density  for  a circu- 
lar coil  probe  for  any  position  and  orientation  of 
the  probe  in  the  dipole  magnetic  field.  This  devel- 
opment allows  the  search  for  Afimaid  to  consider 
“all”  possible  rotations  of  the  three-axis  probe. 
The  extended  search  is  shown  to  leave  unchanged 
the  values  of  ASmax3  that  were  determined  by  the 
earlier  treatment. 


2.  Expression  for  Average  Magnetic  Field 

In  the  derivation  given  below,  it  is  assumed  that 
the  cross  sectional  areas  of  the  wire  in  the  coil 
probes  and  the  opposing  magnetic  fields  produced 
by  currents  induced  in  the  probes  are  negligible. 
We  also  assume  that  the  three  orthogonally  ori- 
ented coils  of  the  three-axis  probe  have  circular 
cross  sections  of  equal  area.  These  assumptions  ei- 
ther can  be  met  in  practice  or  can  be  taken  into 
account  by  a calibration  process. 

The  average  magnetic  flux  density,  5av,  for  a sin- 
gle circular  coil  probe  with  cross  sectional  areay4  is 
given  by 


(2) 

where  (M  is  an  element  of  probe  area,  n is  a unit 
vector  perpendicular  to  A , and  B is  the  magnetic 
flux  density.  In  spherical  coordinates,  the  magnetic 
flux  density  for  a small  current  loop  of  radius  h is 

[4] 


B ~ cosOur  + s'mOuo,  (3) 

where  /xo  is  the  permeability  of  vacuum,  / is  the 
alternating  current,  and  Ur  and  uo  are  unit  vectors 
in  the  directions  of  increasing  r and  6,  respectively. 
The  assumption  is  made  that  h<^r,  and  the  sinu- 


soidal time  dependence  of  the  field  has  been  sup- 
pressed. The  value  of  Bo  is  given  by  the  magnitude 
of  B [Eq.  (3)].  Figure  1 shows  the  spherical  coordi- 
nates r and  6,  a small  current  loop  at  the  origin  of 
the  coordinate  system,  and  a sketch  of  a three-axis 
probe.  The  center  of  the  probe  coincides  with  the 
origin  of  the  prime  coordinate  system  jc',y',  and/'. 
The  coil  probes  are  labelled  PI,  P2,  and  P3,  have 
unit  normal  vectors  /ii,  nj,  and  /13,  respectively,  and 
are  shown  in  Figure  1 (inset)  for  illustrative  pur- 
poses as  being  in  the  directions  of  prime  coordi- 
nates. The  orientation  of  the  magnetic  dipole  with 
respect  to  the  position  of  the  probe  is  characterized 
by  the  angle  6. 


Fig.  1.  Three-axis  magnetic  field  probe  with  its  center  atx=X(„ 
_V  =0,  and  z =Zo.  A small  current  loop  producing  a dipole  mag- 
netic field  is  located  at  the  origin  of  the  unprimed  coordinate 
system.  The  unit  vectors  «i,  /I2,  and  ns  are  normal  to  the  areas  of 
probes  PI,  P2,  and  P3,  respectively.  Changes  in  the  angle  0 cor- 
respond to  varying  the  orientation  of  the  dipole  with  respect  to 
the  probe. 


For  our  purposes,  it  is  convenient  to  express  B in 
terms  of  Cartesian  coordinates.  The  magnetic  flux 
density  is  then  [2] 


B =1 


.3Cxz 


2/--' 
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where  r = ^nd  k are  unit  vectors 

for  the  Cartesian  coordinates,  and  C is  the  constant 
IMJbVl. 
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The  goal  is  to  develop  an  expression  for  5av  at  an 
arbitrary  point  which  can  be  evaluated  for  any  ori- 
entation of  the  coil  probe.  The  value  of  5av3  can 
then  be  found  by  combining  the  rms  values  of  5av 
from  three  orthogonal  directions  according  to  Eq. 
(1).  The  approach  described  below  for  obtaining 
the  desired  expression  for  5av  is  to  transform  the 
problem  into  the  coordinate  system  of  the  coil 
probe.  In  this  coordinate  system,  the  unit  vector 
normal  to  the  plane  of  the  coil  coincides  with  the 
“z-axis”,  B is  expressed  in  terms  of  the  probe  coor- 
dinates, and  the  integration  over  the  area  of  the 
circular  coil  probe  is  carried  out  numerically  in  po- 
lar coordinates. 

We  begin  by  considering,  without  loss  of  general- 
ity, a three-axis  coil  probe  with  its  center  at  x =jro, 
y =0,  and  z=2o,  where  Xo  = rsin0  and  Zo  = rcos6 
(Fig.  1).  We  then  focus  on  coil  probe  PI  and  its 
unit  vector  tti  after  it  is  rotated  through  angles  ai 
and  aj  with  respect  to  the  prime  coordinate  system 
as  shown  in  Fig.  2.  The  unit  vectors  «2  and  m will 
also  change  in  orientation  to  maintain  their  orthog- 
onal relationship,  but  are  not  shown  for  purposes 
of  clarity.  In  the  prime  coordinate  system,  «]  is 
given  by 

/ii  = isinai  cosQ:2+Jsino'isina2  + Arcosai.  (5) 

By  examination, 

/i2  = isin(ai  -l-90°)cosa2+ysin(ai  -I- 90°)sina;2 
-l-itcos(ai-t-90°) 

= i cosafiCOsa2  + j cosaisina2  — Arsinai.  (6) 

The  remaining  unit  vector,  /13,  is  given  by 
n3  = ni  Xn2 

= -isina2+ycosa2,  (7) 

but  it  also  can  be  determined  by  examination,  i.e., 
ai  is  replaced  by  90°  and  02  is  replaced  by  ai  + 90° 
in  Eq.  (5).  For  this  case,  /13  is  constrained  to  lie  in 
thex'-y'  plane.  Later  this  constraint  is  removed. 

The  coordinate  system  of  the  probe  is  reached  by 
the  following  transformations: 

(i)  translation  of  the  origin  to  the  origin  of  the 
prime  coordinates  shown  in  Fig.  1, 

(ii)  rotation  of  the  prime  coordinates  through 
angle  aj  about  the  z'-axis,  yielding  the  double- 
prime coordinates  z",  as  shown  in  Fig.  2,  and 

(iii)  rotation  of  the  double-prime  coordinates 
through  angle  ai  about  the  y"-axis  yielding  the 
triple-prime  coordinates  x'",_y"',z"'  (Fig.  2).  In  the 
triple-prime  coordinate  system,  the  normal  vector, 
til,  given  by  Eq.  (5),  is  along  the  z"'-axis  as  desired. 


Fig.  2.  (a)  Geometry  of  unit  vector  tti  and  coordinates  after 
rotation  of  the  prime  coordinates  through  angle  02  about  the 
z'-axis  and  after  rotation  of  double-prime  coordinates  through 
angle  ai  about  y "-axis. 

Transformation  (i)  is  given  by 

X =xo+x' 

y =y' 

z=Z()-f-z'.  (8) 

The  first  rotation  of  coordinates  (ii)  is  given  by 

x'  =x"cosa2  — y"sina2 
y ' =x''sina2  +>'"cosa2 

Z'=Z",  (9) 

and  the  second  rotation  (iii)  is  given  by 

z"  =z"'cosai  — x'"sinai 
x''  = z"'sinQ:i  +x'"cosai 

(10) 

From  Eqs.  (8)-(10),  we  have 

X =X()  + (z  '"sinoi  +x  '"cosai)cosa:2  -y  "'sina2 
y ={z  "'sinoi  +x  "'cosa'i)sina2  +y  "'cosa2 
z =Z()-l-z"'cosai  — x"'sina:i,  (11) 

which,  when  substituted  into  Eq.  (4),  expresses  B in 
terms  of  the  probe  coordinates.  A simplification  of 
Eq.  (11)  is  had  by  noting  that  the  integration  over 
the  area  of  the  probe  [Eq.  (2)]  occurs  in  thex"'-y"' 
plane,  i.e.,  z"'  = 0.  Contributions  to  Bav3  from  each 
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of  the  coil  probes  is  found  by  using  the  appropriate 
normal  vector  [Eqs.  (5)-(7)]  during  the  integra- 
tions. As  noted  above,  the  integration  is  carried  out 
numerically  (using  a double  Simpson’s  Rule)  in  po- 
lar coordinates,  i.e., 

x"'  = pcosip,  O^p^a,  0^  ijj  <2Tr 

y"'=P  simA, 

dA=dx’"dy"'=p  dp  dip.  (12) 

The  accuracy  of  the  numerical  integrations  was 
checked  by  increasing  the  number  of  divisions  be- 
tween the  limits  of  integration  for  p and  ip.  The 
results  reported  below  were  not  affected  by  further 
refinements  of  the  intervals  used  during  the  inte- 
grations. 

In  the  search  for  A5max3,  it  will  be  necessary  to 
perform  rotations  of  the  three-axis  probe  about  the 
2 '"-axis  or /ii  direction  (see  Search  Protocol  below), 
i.e.,  the  unit  vectors  nj  and  /13  for  probes  P2  and  P3 
are  rotated  about  ni.  This  removes  the  constraint 
noted  earlier  that  lies  in  thex'-y ' plane.  Because 
the  integrand  for  Bav  is  in  terms  of  the  a angles, 
relationships  must  be  found  between  the  angle  of 
rotation  about  the  2 '"-axis,  designated  as  cp,  and 
the  a values  that  appear  in  the  integrands  for  P2 
and  P3.  These  relationships  are  found  by  examin- 
ing the  unit  vectors  for  the  probes  /12  and  /13  as  they 
rotate  about  the  2 "'-axis  or  /ii  direction. 

In  Figs.  3 and  4,  consider  «i  in  a direction  char- 
acterized by  the  angles  am  and  azn,  i.e.. 

Ml  = /sinai()COSa2o+JsinaioSina2o  + ArcosQ'i().  (13) 
From  Eqs.  (6)  and  (7),  we  have 

Hi  = i sin(aio  + 90°)cosa2o  + Jsin(aio  + 90°)sina2o 
-l-itcos(aio  + 90°),  (14) 


and 

Hi  = i sin(90°)cos(a:2o  + 90°)  + j sin(90°)sin(a2o  + 90°) 
+ kcos(W).  (15) 

The  trigonometric  expressions  in  Eqs.  (14)  and 
(15)  are  not  simplified  in  order  to  aid  the  reader  in 
seeing  the  relationships  between  the  three  unit  vec- 
tors. 

Following  a counterclockwise  rotation  of  cp  de- 
grees about  the  2 "'-axis,  the  a2’s  will  increase  in 
value  and  the  ai’s  will  decrease  in  value  in  the  ex- 
pressions for  /12  and  /13.  These  changes  also  occur  in 
the  expression  for  the  magnetic  flux  density  B. 


Fig.  3.  Geometry  of  coordinates  and  unit  vectors  after  unit  vec- 
tor «2  is  rotated  4>  degrees  about  wi  or  r"'-axis.  The  rotation  of 
iti  is  not  shown  for  purposes  of  clarity  (see  Fig.  4). 


Fig.  4.  Geometry  of  coordinates  and  unit  vectors  after  unit  vec- 
tors fi3  and  /12  are  rotated  through  angle  about  «i  or  2 '"-axis. 
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After  a rotation  of  </>  degrees,  a line  along  the 
unit  vector  ni  will  intersect  the  circle  of  rotation  in 
XhQx"'-y"'  plane  at  a point  given  by  (Fig.  3) 

X'"  =CCOS(f) 

y =csin<^ 

z"'  = 0,  (16) 

where  the  radius  for  the  rotation  has  been  arbitrar- 
ily taken  to  be  some  constant  c.  From  Eqs.  (10) 
and  (16),  the  same  point  in  the  double  prime  coor- 
dinate system  is 


z"=  — x"'sinaio=  — ccosc^sinaio 

x"  =x  '"costtio  = c cos</)cosaio 

y"  =y"' =csin(f>.  (17) 

The  increment  to  a2o  for  unit  vector  tij  after  rota- 
tion (f),  §22  (Fig.  3),  can  be  found  from  the  expif^^- 
sion  for  its  tangent,  i.e.. 


X cosaio 


§22  = tan' 


/ tan0  \ 

\ cosaio/’ 


(18) 


Prior  to  the  rotation,  the  angle  with  respect  to  the 
2 '-axis  for  the  unit  vector  tii  is  given  by  am + 90° 
(Fig.  3).  After  the  rotation,  the  corresponding  an- 
gle will  be  §12  + 90°  where  §i2<aio.  The  value  of  §12 
is  found  by  noting  that  after  the  rotation  0,  the  line 
/ from  the  origin  to  the  projection  of  the  circle  of 
rotation  onto  thex''-y"  plane  is  given  by 

l = ^{x"f  + {y")\  (19) 


and  that  the  tangent  for  §12  is  just  |2"//|  (Fig. 3). 
From  Eqs.  (17)  and  (19), 


§12  = tan' 


= tan 


■'( 


cosc^sinaio 


Vcos^</)cos^aio  + sin^cf) 


Thus,  following  a rotation  of  degrees  about  ni, 
a2o  and  am + 90°  will  be  replaced  by  a2o  + §22  and 
§12  + 90°,  respectively,  in  the  expressions  for  /12  and 
B during  the  calculation  of  Bav  for  probe  P2. 

The  a values  for  can  be  determined  with  a 
similar  analysis.  Following  a rotation  of  <j)  degrees 
about  the  2 '"-axis  (Fig.  4),  a line  along  the  unit 
vector  /13  will  intersect  the  circle  of  rotation  in  the 
x"'-y"'  plane  at  the  point 


x'"  = -csin0 
y’"  =ccos<^ 

2'"-0.  (21) 


From  Eqs.  (10)  and  (21),  the  same  point  in  the 
double  prime  coordinate  system  is 


2"  = csin0sinam 
x"  = — csin^cosam 

y"  = ccos0.  (22) 


Following  rotation  </>,  the  angle  a2()  + 90°  for  /13 
[Eq.  (15)]  will  increase  by  an  amount  §23  as  shown 
in  Fig.  4.  The  increment,  §23,  can  be  found  from  the 
expression  for  the  absolute  value  of  its  tangent,  i.e.. 


tan§23  = 


y" 


= tan(/)COSam, 


§23  = tan  '(tan^cosam). 


(23) 


Prior  to  the  rotation,  the  unit  vector  «3  makes  an 
angle  of  90°  with  respect  to  the  2 '-axis  (Fig.  4). 
After  the  rotation,  this  angle  will  decrease  by  an 
amount  §13.  The  value  of  §13  is  found  by  noting  that 
after  the  rotation  cf),  the  line  m from  the  origin  to 
the  projection  of  the  circle  of  rotation  onto  the 
x"-y"  plane  is  given  by 


m=V{x"y  + (y”)\ 


(24) 


and  that  the  tangent  for  §13  is  just  z"lm  (Fig.  4). 
From  Eqs.  (22)  and  (24)  we  have 


tan§i3  = — = 


sin(/)sinam 


Vsin^</)cos^aio  + cos^ 


§13  = tan' 


sinc^sinam 


' \/sin^</)COS^am  + cos^0 


)• 


(25) 


Thus,  following  a rotation  of  (f)  degrees  about  /ii, 
a2()  + 90°  and  90°  are  replaced  by  a2o  + 90°  + §23  and 
90°-§i3,  respectively  in  the  expressions  for  «3  [Eq. 
(15)]  and  B during  calculation  of  5av  for  coil  probe 
P3. 


3.  Search  Protocol 

The  search  for  the  largest  difference  between 
5av3  and  Bo,  A/?max3,  for  a given  distance  r from  the 
dipole  proceeds  as  follows: 

(i)  For  a fixed  distance  r away  from  the  dipole, 
and  with  0 = ai  = a2  = O,  the  three-axis  probe  is  ro- 
tated about  the  2 "'-axis  or  tti  direction  in  2°  steps 
(i.e.,  (})  is  incremented  in  2°  steps).  Then  5av  for 
each  coil  probe  is  evaluated  and  combined  accord- 
ing to  Eq.  (1)  for  each  value  of  (f>  to  obtain  5av3,  5av3 
is  compared  with  Bo,  and  the  largest  difference  is 
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saved.  Because  of  the  symmetry  of  the  problem,  a 
total  rotation  of  90°  is  required  to  cover  all  the 
cases  (with  2°  increments). 

(ii)  The  angle  ai  is  advanced  in  5°  steps  and  the 
above  comparisons  are  repeated  as  the  probe  is  ro- 
tated about  the  z"'-axis  or  tii  direction.  The  maxi- 
mum value  of  ai,  without  duplication  of  results  is 
90°. 

(iii)  For  each  value  of  oi,  az  is  incremented  from 
0°  in  steps  of  5°  and  the  above  comparisons  are 
repeated.  Because  of  symmetry  arguments,  a total 
rotation  of  180°  is  required  to  consider  all  the  cases 
without  duplication. 

(iv)  Following  the  above  calculations,  different 
orientations  of  the  magnetic  dipole  are  considered 
by  changing  the  angle  6 in  15°  steps  and  repeating 
steps  (i)  through  (iii).  The  choices  of  increments 
indicated  above  were  found  to  provide  adequate 
sensitivity  for  determining  ABavs- 

(v)  Steps  (i)  through  (iv)  are  repeated  for  differ- 
ent values  of  r. 

A diagram  schematically  indicating  several  posi- 
tions for  Hi,  and  rotations  about  tti,  as  the  above 
protocol  was  carried  out  for  a fixed  value  of  r is 
shown  in  Fig.  5. 


z' 


Fig.  5.  Trajectories  of  «i  during  search  for  AZJmaxj-  The  unit  vec- 
tors «2  and  It}  are  not  shown  but  maintain  their  orthogonal  rela- 
tionship with  III  throughout  the  search. 


4.  Results  and  Discussion 

As  already  noted,  an  earlier  search  for  ABmaxs  [3] 
was  not  as  comprehensive  as  the  one  described  in 
this  paper.  While  the  ratio  ria  and  0 could  be  varied 
without  restriction,  the  rotations  of  the  three-axis 
probe  were  limited  to  simple  rotations  about  the 
x'-,y'-,  or 2 '-axis.  That  is,  it  was  not  possible  to  con- 
sider differences  between  5av3  and  Bq  for  combina- 
tions of  rotations  about  two  or  three  axes.  This 
problem  has  been  overcome  with  a more  general- 
ized expression  for  Bav  compared  to  the  ones  used  in 
the  earlier  calculation.  What  is  perhaps  surprising, 
however,  is  that  the  ABmaxs  values  obtained  with  the 
more  comprehensive  search  protocol  are  the  same 
as  previously  calculated.  That  is,  ABmaxs  is  negative 
and  occurs  for  all  rja  values  when  6 = 90°,  as  previ- 
ously found,  and  correspond  to  the  ABmaxs  values 
determined  earlier  following  simple  rotations  about 
they '-axis  (referred  to  as  “a  rotations”  in  Ref.  [3]). 
Numerical  values  of  ABniax3  are  provided  in  Table  1 
as  a function  of  rja . 


Table  1.  Values  of  A/?n,ax3  as  a function  of  normalized  distance 
rIa  from  magnetic  dipole 


ria 

A5„„,3(%) 

3 

-19.6 

4 

-10.8 

5 

-6.9 

6 

-4.8 

7 

-3.5 

8 

-2.7 

9 

-2.1 

10 

-1.7 

11 

-1.4 

12 

-1.2 

13 

-1.0 

14 

-0.9 

15 

-0.8 

5.  Conclusions 

The  present  calculations  have  determined  the 
largest  differences  between  the  resultant  magnetic 
field,  Bav3,  and  the  field  value  at  the  center  of  the 
probe  Bo,  assuming  a dipole  magnetic  field.  These 
largest  differences,  designated  ABmax3,  are  reported 
in  Table  1 as  a function  of  normalized  distance, 
rja , from  the  center  of  the  dipole  and  agree  with 
values  previously  found  after  a far  less  comprehen- 
sive search  [3].  The  quantity,  ABmax3,  can  be  re- 
garded as  the  largest  error  due  to  instrumental 
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averaging  effects.  As  noted  earlier,  because  the  rel- 
ative orientations  of  the  dipole  and  three-axis 
probe  are  not  known  for  a given  ria  under  typical 
measurement  conditions,  there  will  be  a range  of 
possible  differences  between  5av3  and  Bq.  Thus,  ide- 
ally, it  would  be  desirable  to  determine  the  distri- 
bution of  differences  between  5av3  and  Bq  and  treat 
the  problem  using  a statistical  approach,  but  that 
has  been  left  to  a future  calculation. 

Because  the  dipole  field  is  a good  approximation 
of  fields  produced  by  many  electrical  appliances, 
the  information  in  Table  1 should  be  taken  into 
account  when  total  uncertainties  are  being  deter- 
mined during  measurements  of  magnetic  fields 
from  appliances.  For  example,  if  the  resultant  mag- 
netic field  is  to  be  measured  at  a distance  r from  an 
appliance  with  a combined  relative  standard  uncer- 
tainty [5]  of  less  than  ± 10%,  magnetic  field  meters 
with  three-axis  probes  having  radii  a such  that 
rIa  ~2)  should  be  considered  unsuitable.  Three-axis 
probes  having  radii  such  that  rja  =5  would  conser- 
vatively be  considered  suitable  if  the  combined  rel- 
ative standard  uncertainty  from  all  other  sources 
(e.g.,  calibration  process,  frequency  response) 
amounted  to  about  3%  or  less,  since  6.9% -t- 
3.0%  = 9.9%,  where  6.9%  is  taken  from  Table  1 for 
ria  =5. 
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